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Abstract- We show that a simple modification of Formulation 1 of Farassat results in a new 
analytic expression that is highly suitable for broadband noise prediction when extensive 
turbulence simulation is available. This result satisfies all the stringent requirements, such as 
permitting the use of the exact geometry and kinematics of the moving body, that we have set 
as our goal in the derivation of useful acoustic formulas for the prediction of rotating blade 
and airframe noise. We also derive a simple analytic expression for the autocorrelation of the 
acoustic pressure that is valid in the near and far fields. Our analysis is based on the time 
integral of the acoustic pressure that can easily be obtained at any resolution for any observer 
time interval and digitally analyzed for broadband noise prediction. We have named this result 
as Formulation 2B of Farassat. One significant consequence of Formulation 2B is the deriva- 
tion of the acoustic velocity potential for the thickness and loading terms of the Ffowcs 
Williams-Hawkings (FW-H) equation. This will greatly enhance the usefulness of the Fast 
Scattering Code (FSC) by providing a high fidelity boundary condition input for scattering 
predictions . 

1- Introduction 


Broadband noise prediction of rotating machinery and airframes is an important problem of 
aeroacoustics. The most successful methods today are semiempirical [1,2] that depend on 
some measured data to find the constants required in the model. We believe that semiempiri- 
cal methods will be useful for the foreseeable future because of their ease of usage and speed 
of execution on a computer. Broadband noise is generated by the random fluctuations of the 
surface pressure as well as velocity fluctuations on the data surface when the governing equa- 
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tion of noise generation is FW - H pds . Here the subscript pels stands for penetrable (permeable 
or porous) data surface. Turbulence in the boundary layer as well as in the medium around the 
moving body, therefore, plays a significant role in generation of the broadband noise. So far, 
the detailed information on the surface pressure and velocity fluctuations has not been avail- 
able. However, we feel that the trend today is toward heavy use of computers for turbulence 
simulation in the future [3- 7]. We, therefore, are going to assume that we have extensive 
turbulence simulation at our disposal and ask ourselves what would be a suitable acoustic 
formulation to use for prediction of broadband noise of rotating machinery and airframes. We 
believe that we can answer this question fully in the present paper. Surprisingly, the answer 
lies in the use of a time domain method. There is also an unexpected byproduct of our analysis 
that has an impact in acoustic scattering computations. 

Time domain methods have rarely been used in broadband noise prediction in the past. We are 
only aware of a few publications by Casper and Farassat [8- 9] on this subject. These authors 
used a time domain result called Formulation IB of Farassat [8]. They then reproduced two 
sets of broadband noise predictions by Amiet in frequency domain for which experimental 
acoustic data were available [10, 11], Comparison of acoustic predictions in time domain with 
both of Amiet’ s predictions and measured data was excellent. Although these studies by 
Casper and Farassat verified the feasibility and suitability of time domain analysis in broad- 
band noise prediction, Formulation IB suffers from some restrictions that do not meet our 
stringent requirements for derivation of a useful acoustic formula for noise prediction. Below 
we present these requirements and the weaknesses of Formulation IB. We will then demon- 
strate that it is possible to derive a new time domain result for broadband noise prediction that 
meets all of our requirements for a good predictive formula. 

To improve broadband noise prediction, a physics based methodology based on FW-H or 
FW - H pds equations is proposed here assuming that the detailed pressure and velocity fluctua- 
tions are known on or near the blade (FW-H eq. [12]) or on the data surface (FW - H pds eq. 
[13]). The Aeroacoustics Branch of NASA Langley Research Center has a long history of 
deriving formulations for rotating blade and airframe noise prediction in the time domain 
(Formulation 1, 1A, etc.). The first author as well as Kenneth S. Brentner (now a Professor of 
Engineering at Penn State University) have been involved in most of these developments. 
Throughout these years, NASA Langley researchers followed a technical philosophy that has 
proven to lead to highly useful acoustic formulations. A set of stringent requirements for 
acceptance of a good predictive acoustic formula underlie the Langley philosophy. These 
requirements are: 

1- There must be no restrictions on the geometry of the noise generator (blades, airframe, 
etc.). This means that results for infinitely thin blades are not acceptable to us. 

2- There must be no restrictions on the kinematics of the noise generator. This means that the 
blade cannot be restricted to a straight line motion or to lie on a helicoidal surface. 

3- The result must be valid in the near and far fields. 

4- One should be able to calculate the noise for an observer that is stationary in the medium or 
in motion with the aircraft. 

We add one more requirement for an acoustic formula to be suitable for statistical analysis 
after all the above requirements are met. We would like the formula to be composed of one 
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term. This can only be appreciated from the authors’ work reported in [14]. From this work, 
one can see that developing statistical results based on Formulation 1A of Farassat [15, 16] 
when the near field terms are retained becomes very complicated even when the body moves 
rectilinearly. The problem can be traced to the fact that we have multiple terms in Formulation 
1 A. We also conclude that Formulation IB does not meet our requirements because it is for an 
infinitely thin airfoil moving rectiliearly with the observer located in the far field. Neverthe- 
less, the assumptions we used in deriving Formulation IB corresponded exactly with the 
assumptions that Amiet used in his analysis in frequency domain [10, 11]. The work by 
Casper and Farassat proved that time domain methods can be used for broadband noise predic- 
tion of rotating blades and airframes. 

2- Derivation of Formulation 2B 


Our governing equation is the FW-H or FW - H pds equation without the quadrupole term. Let 
the moving surface be defined by /( x, t) = 0, / > 0 outside the body. Here the body can be a 
solid surface, e.g., a wing or a blade, or a penetrable (permeable, porous) data surface. Let p' 
be the density perturbation and define p' - c 1 p ' , with c being the speed of sound in the undis- 
turbed medium. The FW-H equation is [12]: 


□ 2 c 2 p' = □- p' = — [p 0 V„ 6(f)] 
at 
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d 
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[p Hi 6(f)], 
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and the FW - H pds is [13] 


d d 

□ 2 p' = — [po u n ] 6(f) - — [Li 6(f)]. 
at oxi 

In Eq. (2), we have used the following two symbols defined by the relations 
P tin 
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(4) 


In these equations, po is the density of the undisturbed medium, v n and p are the normal veloc- 
ity and pressure on the body, respectively. We have used the symbols p and u n as the fluid 
density and velocity with respect to the frame fixed to the undisturbed medium, respectively. 
The component of the unit outward normal to the moving body is 6(f) is the Dirac delta 
function and by is the Kronecker delta. 


Without loss of generality, we work with the solution of FW-H equation retaining the thick- 
ness term involving v„, the local normal velocity of the surface, which in the case of using 
FW - H pds equation, in the form of the term U n of Eq. (3), contributes to the broadband noise 
because this term includes the random fluctuations of fluid particle velocity on the data sur- 
face. Formulation 1 of Farassat is a solution of Eq. (1) for a body in subsonic motion [15, 16]: 
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Here the observer space-time coordinates are (x, t), p = p (y, r) is the local surface pressure at 
the source space-time coordinates (y, r), and r - \x - y| is the distance between the source and 
the observer at the emission time. The symbol M r stands for the Mach number of the source in 
the radiation direction at the emission time. The angle 0 is between the direction r = x - y and 
local normal to the surface. Figure 1 shows the description of some of these variables. The 
subscript ret in Eq. (5) stands for the retarded or emission time which in general is obtained 
from the solution of the equation c(r - t) + \x - y (rj, t)\ - 0 , (x, t) is kept fixed here and 77 is 
the Lagrangian variable of the source on the blade or airframe for which the emission time is 
computed. It is important to recognize that the observer and source variables, x and y, are 
specified in a frame fixed to the undisturbed medium while the variable tj is generally speci- 
fied in a frame fixed to the blade or airframe. 



Figure 1- A moving blade or wing as seen by an observer fixed to the medium at the emission 
time of the source at y. Note that the position of the source in a frame fixed to the moving 
body is described by 77 . 


When a moving body is in uniform rectilinear motion, as is often the case in airframe noise 
problems, the retarded time can be analytically calculated by a geometric construction known 
as the Garrick triangle shown in Figure 2 below. This can result in considerable savings of 
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computation time. 


Radiation process as seen from a frame fixed to the medium. 
TriangleABC is known as the Garrick Triangle. 



A: Source at EmissionTime z 
6: Source at ReceptionTime t 


Figure 2- Computation of emission time and distance using the Garrick triangle. Note that we 
generally know the visual distance r v and the visual angle if/. The radiation distance r can be 
computed by application of cosine law to the triangle ABC from which At = r / c can be com- 
puted to get the emission time t - At. 


We will now write Eq. (5) as follows: 
4 n p'(x , t ) = 
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Here we assume that in the inner integral of the second term on the right, the observer time is 
denoted as t' and that the body started to move at the time t' = 0. We now write the acoustic 
pressure in the following form: 

dU (x, t ) 

P\x, t) = — ^ (7) 

at 


where we have defined a new symbol as follows 
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We will call Eq. (7) Formulation 2B (B for broadband) of Farassat. We propose to use rectan- 
gular integration in time in the second integral of Eq. (8) because for broadband noise calcula- 
tions, the observer time step A t' is a small fraction of the period of the highest frequency of 
interest. The rectangular integration in time that we propose is consistent with the current 
practice of digital signal analysis [17]. We note that the function IT (jt, t) can be calculated by 
any noise prediction code such as WOPWOP [15] or ASSPIN [18] just as the acoustic pres- 
sure is with little computational effort. It is important to note that we do not intend to use Eq. 
(7) for calculation of the acoustic pressure signature which in the case of broadband noise is 
devoid of useful information without further analysis. We propose that the quantity II (jt, t) be 
computed numerically for a given observer time interval at discrete intervals and advanced 
signal processing techniques be applied to these digital data to get the spectral properties of 
the broadband noise. Figure 3 shows the flow diagram of the proposed method for the loading 
term of the FW-H equation. 




BB 

Noise 


Figure 3- The flow diagram of the proposed method of using Formulation 2B shown only for 
the loading noise term of the FW-H equation 


Finally, note that Formulation 2B satisfies all of the requirements of a good noise prediction 
formula that we mentioned in the Introduction. In particular, near field computation of broad- 
band noise is possible from this formula as well as the inclusion of the motion of the blade or 
the airframe. The most striking property of Formulation 2B is that the input data does not have 
to be differentiated in time (as in Formulation 1A) or in space (as in Formulation IB and the 
supersonic Formulation 3 used in ASSPIN). This means that the turbulence simulation data 
can be used without further manipulation and, thus, the accuracy and speed of noise prediction 
are increased significantly. We do not recommend Eq. (7) to be used for discrete frequency 
noise prediction of rotating blades and airframes for which Formulation 1A is recommended 
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[15,16]. 

To find the spectrum of the acoustic pressure signature, we take the Fourier transform of Eq. 
(7) resulting in: 

p'(x,f) = 2niffl(x,f) (9) 

where / now stands for the frequency of the sound and we have defined the Fourier transform 
of n (x, t ) as follows 

r*oo 

n (*,/)= I e- 2nift U(x, t)dt (10) 

kJ — OO 

The fast Fourier transform (FFT) is the most appropriate tool to be used here on the digital 
data IT (x, t). 

We note that the surface terms of FW - H pds equation are similar to those of FW-H equation. 
The fluctuation of fluid velocity by turbulent eddies are included in the U n term of the 
FW - H pds equation. These velocity fluctuations can contribute significantly to the broadband 
noise generation. 

3- Statistical Analysis 


We find the autocorrelation of the acoustic pressure based on Eq. (7). We have, for a suffi- 
ciently large time period T 
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Here R r y (x, r) is the autocorrelation of the acoustic pressure and /?n(x, r) is the autocorrela- 
tion of the digital data n (x, t). We have assumed that we have a stationary random process 
here which is a reasonable assumption in broadband noise prediction. The variable r in Eq. 
(11) is the time displacement of the autocorrelation function and is not related to the emission 
time. It is interesting to note that, in general, the autocorrelation function Rn(x, r) is a smooth 
function for broadband noise so that any good numerical differentiation method can be used to 
find the autocorrelation of the acoustic pressure. 


It is obvious that more statistical analysis can be performed on Formulation 2B. For example, 
the spectral density of the broadband noise can be found by taking Fourier transform of Eq. 
(7) [17], In all these analyses, one always works with the digital data n(x, t ). Ideally, the 
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prediction of the broadband noise as proposed here should correspond to the digital analysis of 
the measured acoustic data, thus, reducing the possibility of errors due to the differences in the 
interpretation of measured and predicted acoustic data. 

4- Further Analysis and Future Plans 


One surprising result following from Formulation 2B is that we have now access to the acous- 
tic velocity potential in the frame fixed to the medium in the near and far fields. Let <p(x,t) be 
the acoustic velocity potential. Then, from the relation 


dll (x, t ) d <p{x, t ) 

P\x, t) = = -p 0 


dt 


dt 


( 12 ) 


we see that 


1 

(p{ x, t) = — n (x, t) (13) 

Po 

This is the acoustic velocity potential for the thickness and loading terms of the FW-H as well 
as FW - H pds equations. To the authors’ knowledge, it is the first time that this result has been 
derived. This is the subject of another paper by the first author, Mark H. Dunn and Ana F. 
Tinetti to be submitted for publication. The acoustic velocity is required to satisfy the bound- 
ary condition in the computation of the scattering of propulsion and airframe noise from air- 
craft wings and fuselage. The Fast Scattering Code (FSC) of NASA Langley developed by 
Dunn and Tinetti [19, 20] can now use this acoustic velocity potential in the near and far fields 
when the incidence noise is computed from any code based on FW-H or FW - H pds . 

The authors are planning to use the results of this paper to reproduce Amiet’s acoustic predic- 
tions of broadband noise in the frequency domain [10, 11], Some excellent measured acoustic 
data by Amiet and coworkers are available in these publications. Following these calculations, 
the authors will study broadband noise from open rotors, helicopter rotors and airframes. 


5- Concluding Remarks 


In this paper we have derived a new formula for the prediction of broadband noise of rotors 
and airframes that satisfies all the stringent requirements we have set for a good noise predic- 
tion formula. We have named it Formulation 2B and it is strictly for broadband noise predic- 
tion when extensive turbulence simulation data are available. The formulation involves the 
calculation of a new quantity n (x, t), the time integral of the acoustic pressure, that can be 
easily calculated by any of our noise prediction codes at NASA such as WOPWOP, ASSPIN 
and ANOPP [21], This quantity can be calculated in any given observer time interval for small 
time increments. Advanced digital signal analysis will then be applied to n (x, t) correspond- 
ing precisely to the analysis of the measured acoustic data to predict the broadband noise. An 
important fact to remember is that no time or space derivatives of the surface pressure or 
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velocity fluctuations are required in Formulation 2B and the turbulence simulation data are 
used as computed. As a matter of fact, an acoustic code can be developed that uses the grid 
structure of turbulence simulation code and thus avoiding the interpolation of turbulence data 
for the acoustic calculation. 
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